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In the framework of the quasipotential method the covariant expression for the 
two-particle vertex operator is obtained. The nuclear structure corrections of orders 
(Za)'^, (Za)^ including recoil effects to gyromagnetic factors of the bound electron 
and muon are calculated. Numerical value of the contribution of order (Za)^ is 
obtained by means of the dipole parameterization for the nuclear charge form factor 
in the range of the nuclear charges Z = 7 32. 
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I. INTRODUCTION 

Experimental and theoretical investigations of the g-factor of a free and bound leptons 
(the electron and muon) are important both for the check of the Standard Model and 
determination more correct values of several fundamental parameters of the theory including 
the lepton-to-proton mass ratio The growth of the interest to theoretical calculations 

on this problem in the last years is connected with new experiments where the hydrogen-like 
ions with moderate values of the nuclear charge Z are used ^ (see experimental data shown 
m Table P). An increase of the experimental accuracy for the electron g-factor measurements 
in such simple atoms generates a need for the calculation of new corrections which can 
become essential for the comparison of the theory and experiment. 

Gyromagnetic factor of the bound lepton can be represented in the form: 

gi = 2 + Agrel + Agrad + Agree + Ag^tr + •••• (l) 

Relativistic corrections Agrfd^ radiative corrections Agrad, recoil corrections Agree were in- 
vestigated in Refs.|3, H, IE M H S HO, 11, E^; with an accuracy up to fourth order 
terms. Certain contributions to the gi of higher order over a related for the most part to 
the vacuum polarization and the lepton self energy were studied in Refs.|il El Q. An 



agreement between theoretical and experimental results in the problem of the gyromagnetic 
factor of the bound electron is observed to the present. For hydrogenic ions with sufficiently 
high values Z the nuclear structure effects Agstr acquire important role. They can be taken 



Talk presented on poster session of the Bogolubov Conference " Problems of Theoretical and Mathematical 
Physics", Moscow-Dubna, September 2-6, 2004 



2 



TABLE I: Experimental data and theoretical values on the g-factors of the bound electron. 



Ratio 
of the g-factors 


Experiment 


Ref. 


Theory 
[1, 3] 




i — i/./uy(^ioj X iu 


1 ioj 


i — LI. ( UOo X iU 








1 — 17.7125 X 10 








1 _ 17 501 X 10-6 


(H) 

g.-{D) 


1 + 7.22(3) X 10-9 


\19] 


1 + 7.247 X 10-9 


g.-(T) 


1 + 10.7(1.5) X 10-9 


[20] 


1 + 10.7 X 10-9 


ge-rHe+) 


1 - 70.87(30) X 10-6 


[21] 


1 - 70.87 X 10-6 


g,-(^^C^+) 


1 + 520.798(2) X 10"^ 


[22] 


1 + 520.796 X 10-6 


ge-('"0' + ) 

2 


1 + 23.514(2) X 10-6 


[23] 


1 + 23.511 X 10-6 



into account by electromagnetic form factors of minimal multipolarity describing the distri- 
butions of the electric charge and magnetic moment. Main contribution of order {Za)"^ to 
the Aggtr is determined by the nuclear charge radius Rn which is the differential parameter 
of the electric charge distribution. It was originally obtained on the basis of nonrelativis- 
tic approach in Ref. and by means of the Dirac equation in Ref. 0| where relativistic 
correction of order {Za)^m{B^^ to the ^Qstr was found also. The quasipotential method 
for the problem of magnetic moment of hydrogenic atom was formulated in Refs.|3, 13 • 
this approach the two-particle vertex function describing the interaction of the bound state 
with external electromagnetic field was introduced. Its expansion by the perturbative series 
for hydrogenic systems allows to calculate radiative, recoil corrections of different orders. 
In the case of the constituent particles with arbitrary spin the construction of the vertex 
operator in two main orders by the perturbative series was carried out in Ref.^3|- The aim 
of present work consists in the calculation of the nuclear structure effects of orders {ZaY 
and {Za)^ including recoil effects in the g-factors of the bound electron and muon on the 
basis of the quasipotential method. Considering that in recent experiments with hydrogenic 
ions an accuracy comprises several ppb [3, IT3, H^, 21, [i^, H^J it is evident that the role 



of the nuclear structure corrections in the lepton g-factor for moderate and high Z will be 
enhance. 



II. GENERAL FORMALISM 

Magnetic moment of the two-particle bound state can be written as follows 

d 



M = -'- 
2 



(9A 



A = Q - P, (2) 



where the matrix element of electromagnetic current operator J between the bound states 
with total momenta Q and P respectively can be expressed in terms of the wave functions of 
the bound system \l/„.Q(q), \l/„^p(p) and generalized vertex function represented in Fig.l 
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Pi 



N > — #^ ^ ^ 

(a) (b) 

FIG. 1: Generalized two-particle vertex function F^: the diagram (a) represents F^''^ with free 
noninteracting particles, the diagrams (b) determine the nucleus (A'^) structure corrections of order 
(Za)^ in the g-factor of the bound lepton (l). Bold line denotes negative energy part of the lepton 
propagator. Large bold circle denotes the interaction vertex of the nucleus with electromagnetic 
field proportional to the factor (Fi(A;^) — l). 

in the form: 

< ^„,q|J^(0)|^„,p >= J ^^„,Q(q)r^(p,q)^^„,p(p), (3) 

The two-particle vertex function is determined in terms of the five-point Green-like 
function 

R, =< o\Mt,^i)Mt,^2)Mo)Mr,yi)Mr,y2)\o >, (4) 

projected on the positive energy states 0, 

F^ = = uiU2R^u,U2, (5) 

where G is the two-particle Green function and Ui, U2 are the wave functions of a free 
particles. For loosely bound system in quantum electrodynamics all introduced quantities 
F, R can be represented by the perturbative series: 

F = fW + F(i) + F(2) + = i?o + ^1 + ^2 + = Go ' - 1^1 - .... (6) 

Substituting Eq.(6) in Eq.(5) we obtain the following relations for generalized vertex func- 
tion: 

FW = G,^R^G^\ F(i) = Go ii?iGo ^ - FiGoF^") - F^GqI^i, (7) 

F(2) = Go 'i?2Go ' - l^2Gor(°) - F(°)Go\^2 - ^iGor(i) - F^^^GoV^i - 2ViGoV^''^GoVu 

where Go is the Green function of two noninteracting particles, Vi is the quasipotential of 
the one-photon interaction, V2 is the quasipotential of the two-photon interaction. 

The transformation law for the wave function \l/„^p(p) of two bound particles with spins 
Si, S2 from the rest frame to the reference frame moving with momentum P was obtained 
in Ref. 3 in the form: 



vl/„,p(p) = DnRZ)DnRt)'^nM. (8) 
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where D'^lR) is the rotation matrix, is the Wigner rotation and Lp is the Lorentz boost 
from the rest frame to the reference frame moving with momentum P. The rotation matrix 
can be expressed in terms of the Lorentz transformation matrices as follows: 

D'iRt) = S-\pi,2)S{P)S{p). (9) 

The quasipotential bound state wave function \E'o(p) in the rest frame of the composite 
system satisfies the quasipotential equation |25l |: 

Co^^ = - ^1 ^o(p) = / V(p,q,M)vl/o(q)-^, (10) 

where fiR is relativistic reduced mass: 

_ EiE2 _ - {ml - mlf _ - ml, + ml^ 
^« - ^ - UP ' - 2M ' ^^^^ 

M = El + E2 is the mass of the bound state, 

In the nonrelativistic limit the equation (10) reduces to the Shroedinger equation with the 
Coulomb potential. The main contribution to the vertex function is equal 

r(o)A.,p.(p^q) = «t(qi)7M<(Pi)^2(q2)<(P2)(27r)35(p2 -q2)5'^^ A,a,p,^ = ±1, (13) 

where the Dirac spinors are 



jiL x^ .^4 (14, 

^ V e(p) +m J 



Relativistic four momenta of the particles in the initial and final states are defined as follows: 

3 p 
Pi,2 = ei^2{p)v ±J2n^'\v)p\ v = —, (15) 
1=1 

1=1 

and n^*^ are three four vectors defined by 

n«(t;) = S'^ + -^t;vj . (16) 

Using relations (8)- (9) we can transform the matrix element (3) to the form containing the 
product of two traces over the spinor indices of both particles. For this aim the following 
relations are useful: 

S^f,{A)4{p) = f: <(Ap)D:,«^), (17) 
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a= — .s 

Substituting expressions (8), (9), (13) in Eq.(3) and using relations (17) we obtain: 
<^.d^°^I^.P>-/^^0(q)..(0)- + 



(2v 

{pi + mi] 



2ei{p) + mi) 

{P2 



=Mi(0)m2(0) 



^2ei(g)(ei(g)+mi) 
(ga + "^2' 



S-\LQh,S{Lp) X (18) 



xS{L 



m2 



2e2{p){e2{p) + m2] 



S-\Lq)x 
2e2{q) + m2) 

M2(0)5(p2-q2)^o(p), 



where 



Pi = (ei(p),P),P2 = (e2(p),-p),gi = (ei(g),q),g2 = (e2(g),-q). 
Introducing the one-particle projectors on the states with definite spins si^2 in the rest frame 

:i+7°) (1+7'M 



7ri,2 



^1,2(0)^1,2(0)] 



(19) 



we can represent basic matrix element (18) in the form: 

dpdq 



< V^n,Q| ji°HO)|^n,P >= / ^^0(q)^0(p)5(P2 " ^2) X (20) 



xTr 



xTr 



{qi + mi] 



^2ei(g)(ei(g)+mi) 
(g2 + m2) 



S-\Lq)^,S{Lp] 



-.S-\LQ)SiLp) 



(pi + mi) (1 +70)(1 +7551) 



^2eiip){ei{p) + mi) 4 

(p2 + m2) (l + 70)(l + 75s2) 



Y2e2(g)(e2(g) + m2) ' ^' ' ' •^2e2(p)(e2(p) + ^2) 

In the next orders of the perturbative series the matrix element of electromagnetic current 
keeps the general structure (20) but the form of the vertex function (5) will be changed. 
The expression (20) is useful tool for the calculation of the different order corrections to the 
lepton gyromagnetic factors by means of the computer systems of analytical calculations. 
The system Form [2^ is employed in the calculation of the correction Agstr- 



III. NUCLEAR STRUCTURE CORRECTIONS OF ORDER (Za)^ 

In the framework of the quasipotential method formulated in the section 2 there are two 
sources of the nuclear structure contributions of order {Za)'^ to the gyromagnetic factors 
of the bound particles. The first contribution is represented on the diagrams Fig.l (b). 
The interaction vertex of the second particle (nucleus) with electromagnetic field (large bold 
circle) is determined by two form factors Fi^2- 

Let us consider the contributions of the diagrams Fig.l (b). The trace over the lepton 
spinor indices for the first diagram Fig.l (b) entering in Eq. (20) takes the form: 



Til 



64mf 



Tr (1 + 7°)(1 + ^'si){qi + ^1)7° \v' + 1 - {1 - 7°) 



(21) 



(c) 

FIG. 2: Generalized two-particle vertex funetion F^: the nuclear structure corrections of order 
(Za)^. Bold line denotes the negative energy part of the lepton propagator. 



xA-(p;)7V + 1 - (1 - 7°)] 7°(Pi + 



A-(p'i) 



2mi 



7 



mi(/-l)+P + A(l-^) 



,pi = Pi + A, p=(0,p). (22) 



Similar contribution of the second diagram Fig. 1 (b) is equal 

= ^Trl{l + 7°)(1 + 7'h){qi + m^h' [v' + i-(i- ^' 



Y X 



(23) 



xA-(q'j7V + 1 - (1 - 7°)1 7°(Pi + ^i) , 



A-(q'i) 



2mi 



7 



m,(/-l) + g-A(l-^) 



,q; = qi-A, g=(0,q). (24) 
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The trace over the spinor indices of the second particle (nucleus) which is equal for both 
diagrams Fig. 1 (b) has the following structure: 



21 



Tr (1 + 7°)(1 + 7's2)(g; + m2)7° f)' + 1 - (1 - 7° 



7°x 



(25) 



2m' 



v + l-{l- 7°)1 (^2 + ms) 



After the contraction of the expressions (21) and (23) with (25) over the Lorentz index u 
{fi i is the vector index) we obtain that the sum of the trace products can be written in 
the form: 



[(T^ X A]Fi 1 



mi 



nil + 1^2 



[0-2 X A] (Fi + 2F2) 



nil 



(26) 



nil + ^2 

Substituting sequentially Eq.(26) in Eq.(20) and then in Eqs.(3) and (2) and calculating the 
roiA we keep only the contributions to the lepton gyromagnetic factor (the particle 1). The 
nuclear structure correction takes the form (si — cri >): 



Agf = -mfRj^iZay 1-4 



mi 
ni2 



where the following expansion of the Dirac form factor Fi at small is used: 

k2 



Fi(k' 



1 



6 



-R 



N- 



(27) 



(28) 



The expression (27) contains not only the correction of order {Za)* but also the recoil 
effects of the first order over the ratio nii/ni2. The same order correction as (27) appears 
in the second order of the perturbation theory. The contribution of the second order of the 
perturbative series to the matrix element (3) has the form: 



2 < ^„|J(0)|5^„ >= 2 < ^„|JGAF|^„ >= 



< ^n|J(0)|^^ >< ^^lAl^l^n > 

rp rp ' 



(29) 



where are the Coulomb wave functions of Eq.(lO). The reduced Coulomb Green function 
(RCGF) G doesn't contain the free two-particle Green function because its contribution 
cancels the term which appears from the half of last iteration contribution T^"^^ in Eq.(7). 
The nuclear structure effects in Eq.(29) are determined by the terms of the quasipotential 
AV ~ R% as it follows from Eq. (28). To calculate the matrix element in Eq.(29) we need to 
know the form of relativistic correction of order p^/m^ including recoil effects in the vertex 
operator which was found in Refs.0, HI- Accounting the form of the IS state RCGF in 
the coordinate representation ji^] 



Gis{r,0) 



X 



2x(lna; + C) + — 5x 



(30) 



(C=0.5772156649... is the Euler constant) we can express necessary nuclear structure cor- 
rection in the g-factor of the bound lepton in the second order of the perturbative series as 
follows: 



Ag. 



str 
2 



--RljTiZa / ^isir 



3mf 



1 + ^) 
2m2y 



dvGisir, 0)1^, 



(31) 



8 



TABLE II: Nuclear structure corrections in the g-factor of the bound electron and muon of orders 
(Za)^ and (Za)^. 



z 


RNifm) 


6^s[xl0-9] 
\13] 


electron IS.g'^^'' 


muon Agf**'^' 


(Za)4[xl0-9] 


(Za)5[xl0-i2] 


(Za)4[xlO-4] 


(Za)5[xl0-6] 


7 


2.54 


0.79 


0.79 


0.13 


0.33 


-0.17 


8 


2.737 


1.56 


1.56 


0.26 


0.65 


-0.42 


9 


2.90 


2.80 


2.80 


0.53 


1.17 


-0.89 


10 


2.992 


4.54 


4.54 


0.95 


1.90 


-1.66 


12 


3.08 


9.98 


9.97 


2.49 


4.19 


-4.54 


15 


3.191 


26.14 


26.13 


8.10 


11.03 


-15.50 


18 


3.423 


62.37 


62.37 


22.80 


26.39 


-47.88 


24 


3.643 


223.29 


223.28 


117.13 


94.68 


-244.35 


32 


4.088 


888.63 


888.60 


550.02 


377.57 


-1461.07 



Summing two corrections (27) and (31) we obtain total contribution of the nuclear structure 
effects of order {ZaY with the account recoil correction: 

A^fr^.). = iRWriZaf (l - ^) . (32) 

Numerical values (32) for several hydrogenic ions are presented in Table IHl 



IV. SPIN DEPENDENT NUCLEAR STRUCTURE CORRECTIONS OF ORDER 

{ZaY 

In general case, the matrix element of the electromagnetic current for the particle of arbi- 
trary spin S is determined by means of (25" + 1) form factors (charge, magnetic, quadruple, 
et al.). When the magnetic moment of simple atomic systems is studied it may be possible 
to take into account the form factors of the minimal multipolarity describing the distribu- 
tions of the electric charge and magnetic moment. The one-particle matrix element of 
the electromagnetic current operator between states with momenta p and q can be written 
as follows: 

J, = f/(p) jr^F/^ + ^S,.A;^Ff } f/(q). (33) 

The wave function U{p) of a particle with arbitrary spin entering in Eq.(33) can be presented 
in the form (see, for instance |i2S. .2S]): 



r/- 



U={l] = \Z\ti 1' P + ^ = ^S, (34) 



9 



where the spin-tensors ^, t] are symmetrical in upper and lower indices. For the particle 
of half integer spin p=S+l/2, q=S-l/2. In the case of integer s pin p=q=S. The Lorentz 
transformation of the spinors ^ and rj can be written in the form 2^ 30l |: 



C = exp 



^0, V = exp 



S5\ 



Co 



(35) 



where the direction of the vector d coincides with the velocity of the particle, thS = v. The 
generator of the Lorentz transformation S is equal: 



p 

i=l 



p+q 

E ' 

i=p+l 



and cTj acts on the ith index of the spinor follows: 



(36) 



(37) 



In the standard representation, which is introduced in analogy with the spin 1/2, the free 
particle wave function (34) can be written with the accuracy (f/c)^ in the form: 



1 + 



(Sp)2 



8m2 



(38) 



The components of the matrix S^,^ in (33) are the generators of the boosts and rotations 



-"aO 



S„ 

-s. 



-2ie. 



mnk 



Sk 
Sk 



25 

2^ 

1=1 



(39) 



Considering the contribution of the diagrams Fig.l (b) to the matrix element (20) let us 
point out that with desired accuracy the bispinor contraction over line of the second particle 
(nucleus) 

1 



f/2(q2)i32f/2(P2)i^i(A;' 



1 



6 



2 d2 

N 



-k'R 



(40) 



(where the matrix B2 is the generalization for the matrix /32 used in the case of the spin 1/2 
particles) doesn't contain any other additional factors connected with the nuclear spin. So, 
the form of the correction (27) remains unchanged with the account the first order recoil 
effect. But the correction of the second order over the lepton-to-proton mass ratio contains 
the spin dependent terms. Indeed, using Eqs.(33), (38)- (39) we have the following bispinor 
contraction for both particles (see the diagrams Fig.l (b) where the nucleus interacts with 
external electromagnetic field): 



62 

2m2 



Ki(qi)a>i(pi)Fi(P)f/2(q2)S2 



J^^]- (1 - B2) A2 + (1 - B2) A 



U{P2) 



Ze 



mi 



i[ai X ^]Ks,^\^R^ 



2m2 mi + m2 



6 



Ni 



where 



S 
S 



B 



I 
-/ 



(41) 
(42) 
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TABLE III: Spin dependent nuclear structure corrections in the g-factor of the bound electron 
and muon of order (Za)^: Agf = lR%ml{Za)^^ [QKs^ + Z + 29]. 



z 


RNifm) 


.52 


electron [xlO 


muon [xlO ^] 


7 


2.54 


1/2 


0.12 


0.02 






1 


0.13 


0.02 


8 


2.737 





0.17 


0.03 






5/2 


0.27 


0.05 


11 


2.94 


3/2 


0.53 


0.10 






3 


0.63 


0.12 






4 


0.71 


0.13 


13 


3.035 


5/2 


0.94 


0.17 






3 


0.97 


0.18 


24 


3.643 


7/2 


5.28 


0.97 



^ < Eg > I S2 is integer nuclear spin 

3 I ^^3^ ' ■^s is half integer nuclear spin 

The following commutation relations are useful: 

[Si, Sj] = AieijkSk, pi, Sj] = ieijkTjk- (44) 

Other part of the contribution of order {Za)^m1/m2 is determined by relativistic effects in 
the vertex operator F*^'') which were obtained in Ref. [l^ on the basis of relations (34)-(39) 
in the form: 



(0) _ ei J p2 / mi \ p2 Z\ Z-1 



KsJl + - 



3/ 6 



(45) 



Calculating the contribution of the second order of the perturbative series by means of 
Eqs.(29) and (44) and taking into account the relation (41) we obtain the following total 
value of the nuclear structure correction: 

^9t!Lai = lrnlR%iZar {l " ^ + £| + ^ + 29]| . (46) 

It is helpful to remark that two terms in Eqs.(41) and (45) proportional to Kg^ ■ Z will 
be mutually cancelled in the correction (46). The expression (46) is the generalization of 
Eq.(32) for the case of the nuclear spin dependent contribution which contains also the 
squared ratio of the particle masses mi/m2- Numerical values of the third addendum in 
Eq.(46) for different Z and S2 are compared with the correction of order {Za)^mlR% 
and represented in Table UTTl 



V. NUCLEAR STRUCTURE CORRECTIONS OF ORDER (Za)^ 



The nuclear structure corrections of order (Za)^ in the g-factor of the bound lepton are 
specified by the vertex function F^^^ in Eq.(7) and shown in Fig.2. There is need to take 
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into account that the corrections of such order contain also additional small factor mi /A 
(A is the form factor parameter) which lead to the decrease this contribution as compared 
with Eq.(32). Nevertheless just as for the g-factor of the bound electron in hydrogenic ion 
with sufficiently high so for gyromagnetic factor of the bound muon the value of this 
correction can increase essentially and become important for modern experiments. 

The two-photon exchange diagrams relevant to the T^"^^ can be divided into two parts. 
In the Feynman amplitudes of the first part the interaction of the lepton with external field 
occures between two exchange interactions. We consider these two-photon exchanges in 
the limit of static nucleus neglecting recoil corrections when the nucleus interacts with the 
photons via the Dirac form factor Fi. The interaction of the lepton with two exchanged 
Coulomb photons is plotted in Fig. 2 (a). The bispinor contraction over the lepton line with 
consideration for the interaction terms (7) incorporates negative energy parts of the lepton 
propagators shown in Fig. 2 (a) by bold line. The trace over the spinor indices of the first 
particle entering in Eq. (20) can be written as follows: 

Tu = ^Tr{{l + 7°)(1 + l'h){qi + mi)[v' + 1 - (1 - 7°)]7° X (47) 

X ei(k + A) - A;° + 7°(^ + A-mi)]7V [ei(k) - A;° + 7°(A; - mi)] [{)+l-(l-7°)](pi+mi)} 

f» 4mi(ei(A;) - mi)[cri x A]. 
Then we can write the contribution of the diagram Fig. 2 (a) to the g-factor of the bound 



lepton in the form (ei(A;) = yk'^ + ml): 

(1) 29m,fi^{Za)' p dfc (Ff - l)(6i(fc) - 

327r io k^ {ei{k) + mi)^ei{k)^ ' ^ ^ 

The scale of the integration momenta in Eq.(48) is of order of the lepton mass. In spite of 
the infrared finiteness of the appeared expression in the leading order over mi/m2 we make 
the subtraction of the point-like nuclear contribution in Eq.(48). To obtain correct energy 
spectrum in the case of point-like nucleus it is necessary to consider the bound state effects 
and keep the relative motion momenta of the particles p, q (see the discussion, for example, 
in Ref. [3l|). Another part of the two-photon exchange diagrams which lead to the nuclear 
structure corrections of order (Za)^ is presented in Fig. 2 (b), (c). The trace calculations 
over the spinor indices and the contraction over the Lorentz indices can be performed for 
these diagrams in the same manner as for previous contributions. To illustrate the general 
structure of intermediate expressions in this case we write the trace product regarding to 
the first diagram in Fig. 2 (b): 

rr{(l+7°)(l+7^si)(gi+mi)({)' + l)7^A-(-A)7V(Pi + ^+^i)7"(^ + l)(Pi+mi)}x (49) 

xTr{(l+7°)(l+7^S2)(g2+m2)(^'+l)7^(rn27°-^A-A;+m2)7"(^)+l)(p2+m2)}Fi(fc)Fi(fc+A). 

Total contribution of all diagrams in Fig. 2 (b), (c) in the trace product of Eq.(20) is the 
following: 

f ^0 32 
m2F^\ - — k^^[(Ti X A] - 32mX[cTi x k] - —k^[ai x A] + (50) 
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14 I 
+8k^[a-i X k] - —k'^kllcTi x A] + IQk'^kHai xk]\. 

But contrary to earlier considered expressions for the matrix elements of electromagnetic 
current (2) which are linear with respect to [cri x A] these diagrams involve also the integral 
of the form / d'^k[cri x k]F(fc, A). Its calculation includes the shift of the integration variable 
k ^ k — Y and the calculation of the rot a from the integration function. After rotating the 
k^ contour of the Feynman loop integration we integrate over four dimensional Euclidean 
space using the relation 



d'^k = An 



k dk I sm~ 





k^ = kcosi 



(51) 



The integration over the angle variable (p can be done analytically. After that we can present 
the sum of the contributions shown in Fig.2 (b), (c) to the lepton g-factor in the form: 



(2) 

;,(Za)5 



3 



TT 



dk 



X 



(52) 



X 



TT 



25/3 + I6/1 + Wk'^h + 31/2 + 8k'' 



d 



Vl — x'^dx 



dnii 
UJ2 - 



16Fi'(0)- 



' u f 

mim2k J 



-1 (P + 4mfx2)(A;2 + Amlx"^) 



k"^ + Amf, 



TV 



X^y/1 — x'^dx 



Ak-^^ml — ml) ' * 
k{ml — 7712) + 777,2^71 — mluJ2 



{k^ + Amlx'^){k'^ + Am^x'^^ 



16k{m2 



mDmlml 



(53) 



(54) 



^Vl — x'^dx 



2mlm\{m\ — m^) + k'^{m\ — mg) + A;(m|ci;i — m\uj2) 



TT J-1 (/c^ + 4mfx2)(A;^ + Am^x"^) 



h = - 



QAmfm2{ml — ml 



;i - xy/^d 



X 



ml)^ 



TT J-1 (A;2 + 4mix2)2(/c^ + Am^x"^) 
-k'^ + k'^{Qm\ — 2m\ + UJ1UJ2) + 4(2m^ — Qm\m\ + m\uJiUJ2) 



(55) 
(56) 



IQk^uJiiml 

The contribution of the iteration term in the vertex operator F*^^) (7) was taken into account 
in Eq.(52). We make also necessary subtraction of the point-like nuclear contribution in 
Eq.(52) because in the low momentum region there exists logarithmic infrared divergence. 
To obtain numerical results on the basis of the expressions (48) and (52) the dipole param- 
eterization for the Dirac form factor Fi is used: 



Fi{k' 



1 



A 



12 



R 



N 



(57) 



The total contributions of order {Za)^ to the Agstr for the electron and muon presented in 
Table ITTl have different signs because the corrections (48) and (52) have opposite signs and 
their numerical values differ considerably. 
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VI. CONCLUSION 



The calculation of the nuclear structure corrections (the nuclear size corrections) of orders 
{ZaY and {Za)^ in the electron and muon g-factors is performed in this work on the basis 
of the quasipotential method. Numerical results of obtained corrections (32), (46), (48), 
(52) are presented in Tables |n] and IIII [ We kept the same values of the nuclear charges 
and nuclear charge radii as in Ref.jjj|. It was demonstrated that considered corrections 
can be successfully derived in the framework of the quasipotential method used earlier in 
Refs.j3, 1^ for the study the radiative and recoil effects of orders 0{a^), O(a^). The matrix 
element of electromagnetic current operator is obtained in new form (20) which is convenient 
for the calculation of different quantum electrodynamical corrections of high order over a 
and mi/m2 by means of such systems of analytical calculations as the Form, Reduce etc. 
Obtained results for the corrections of order {ZaY coincide with the calculations in Refs. 
jisi [l6|. New recoil contributions {Za)^mi/m2, {Za)^{mi/m2)^ and the correction of order 
(Za)^ in the g-factor of the bound lepton are numerically less important to the present than 
the leading order contribution. Their role can be significantly increased with the growth 
of experimental accuracy for the hydrogenic ions with sufficiently high Z. All obtained 
analytical relations and numerical results for the correction Ag^^^ in the IS* state can be 
extended to other S states of simple atoms. It should be pointed out that for muonic 
hydrogenic atoms the role of the nuclear structure corrections of orders {Za)^ and {Za)^ 
including recoil effects in the muon g-factor increases by some orders of the magnitude 
as compared with the g-factor of the bound electron. The theoretical expression for the 
gyromagnetic factor of the bound muon in the hydrogen-like atom accounting relativistic, 
radiative, recoil and structure corrections of the fourth order has the form 



9,iH) = 
12 



2{1- 



{Zaf 



2a{Za) 



7r 



2 


1 - 


3 rrif, ■ 






2 misf. 




"2 




'2Zam 






In ( 




9 




- rrie 





'l + Z){Za) 



m„ 



TUN 



2 ajZaf 



1 m,. 



3 



27 



+ 



+ 0{mf,/mN) 



- (58) 



Zaf). 



2{Zam^ 

In the case of the muonic hydrogen numerical value of the nuclear structure correction of 
order [Za)^ is equal 1.11 x 10~^ whereas the total value of the corrections to the gyro- 
magnetic factor QuiH) following from Eq.(58) is (1 — 15490 x 10~^). So, an enlargement 
of experimental region to the g-factor of the bound muon in the muonic hydrogenic atoms 
which is studied in this work, and an enhancement of experimental accuracy can lead to 
more precise determination of several fundamental physical constants such as the proton 
charge radius and the muon-to-ptoton mass ratio. 
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